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FIG. 1. A schematic illustration of a chaotic billiard with 
a rough boundary threaded by AB flux <f). An example of an 
unstable periodic orbit is shown. 



d(E)=d°(E)+d osc (E,< 



(4) 



The Weyl term depends on the classical phase volume 
and is independent of 0. (The higher ft corrections to 
d° give the standard bulk Landau magnetism£3) d osc (E) 
can be written as a sum over classical periodic orbits s, 



d° sc (E,4,) = 



-y 



det(M s - 1)\ 1/2 



X cos 



— - + 2tTUJ s — 
ft 00 



7T 
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where r s , M s , S* s , w s , and a s are the period, the mon- 
odromy matrix, the action integral, the winding num- 
ber and Maslov index, respectively. Following Richtcr, 
Ullmo, and Jalabert,Ej we define a mean chemical poten- 
tial ft by the condition of accommodating TV electrons to 
the mean number of states, 



N = N(p) = 7V(/Z). 



Here 



7V( M ) 



/ dEd(E)f(E-n), 

Jo 



(6) 



(7) 



where f(E — fi) is the Fermi-Dirac distribution func- 
tion. N(p) is obtained in eq (Q) by replacing d(E) and 
// by d°(E) and /Z. The Helmholtz free energy in terms 
of grand canonical quantities is obtained by expanding 
eq. (0) to second order in [i — p, as 



F(N) w F° + AF (1) + AF< 2) , 



(8) 



with 



F° 


= p,N + n(ji), 
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where L s = vpr s is the length of trajectory s, and d(/i) = 
mA/2Trh 2 (m is the mass of electron and A is the area of 
the billiard). In eqs. @ and @, 



i?r(L s ) = 



L s /L c 



sinh(L s /L c )' 



(12) 



with i c = hvFf3/ir, yields an exponential suppression for 
the contribution of trajectories with L > L c . 

The leading-order contribution to F is given by F° + 
AF^\ and yields the persistent current of single chaotic 
billiard calculated in the grand canonical ensemble as 






Rt{L s )uj s 



L s |det(M s -/)| 1/2 
x sin ( kpL s + 2iroj s - a s 



(13) 



where Up = hkp/m is the Fermi velocity. Since the sys- 
tem is time reversal invariant at zero magnetic flux, each 
periodic orbit is associated with a time-reversed part- 
ner having exactly the same characteristics except for an 
opposite winding number. Grouping the time-reversal 
trajectories t (u>t > 0) in eq. (|13|), we have 
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(14) 



The persistent current for single billiards exhibits a 0o 
periodicity and changes from diamagnetic to paramag- 
netic by changing kp. 

In order to characterize the typical value of the persis- 
tent current for single chaotic billiard, we define 



I typ = V /(K 1 )(0) 2 ) 



(15) 



0/00 = 1/4 



where the average (• • •) denotes an averaging over kp or 
some external parameter associated with the roughness 
of the boundary. Upon averaging, the off-diagonal terms 
vanish due to the widely fluctuating phases, while the di- 
agonal term survives the averaging. Therefore we obtain 



/M(0) 2 ) = 



y/2, 
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To evaluate the sum of t, we order the orbits by their 
length L t and use pthe Hannay and Ozorio de Almeida 
(H-OdA) sum rule,ll3 i.e., 



E 



det(M t -I)| 



ff±PU 
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Here Pl{u) denotes the classical distribution of u that 
an orbit is of length i^jand is approximately given by 
Gaussian distribution,! 



Iff 



Pl(lo) 



1 



: exp 



fin^T) V 2 ^ L ), 



(18) 



with a variance increasing linearly with L, i.e. 

L 



u*(L) 



V 



(19) 



where a is a system dependent constant, and Lq is the 
length for the shortest periodic orbit. 

In the following, we calculate P vp in high-temperature 
regimes, T > T^EHl where T t h = hvpn/ksL* is the 
threshold temperature. (L* is the characteristic length 
for which periodic orbits-can be taken as uniformly dis- 
tributed in phase space.Eil) In this regime, the thermal 
fluctuation washes out structures on the energy scale 
of the mean energy level spacing A, and only shor±_or- 
bits contribute effectively to the persistent current .L 2 ] In 
this case one may exploit the asymptotic form R<r(L) sa 
2(L/L c )cxp{-L/L c ). Substituting eqs. @ and @ 
into eq. (pf), and using the Poisson summation formula, 
lead to 
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(20) 



where ui c — ^JaL c j ' Lq. Therefore, the typical value of 
persistent current is given by 



pyp — 



2^2 



evF 



,-l/a 



: c (1 - e -2M) vi + e- 2 ^ 
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(21) 



The typical persistent current [eq. (pl|)] scales as I* w ~ 
fcp^, which is in agreement with th&xesult of von Oppen 
and Riedel's dimenptwaal analysis Hi In the case of in- 
tegrable systems J130E3 on the other hand, it scales as 
pvv ^j (kpt) . This mean that the persistent current 
of integrable billiards is much larger than that of chaotic 

As in the case of diffusivetJ~Q and integrableEj systems, 
/t 1 ) gives a vanishing contribution to the grand canon- 
ical persistent current of an ensemble of billiards with 
different sizes or electron fillings as soon as the disper- 
sion in kpt is order of 27r. Therefore we need to go to the 
canonical term AF^ 2 ' in the free-energy expansion. Us- 
ing eqs. |y) and (|ll|), for the canonical persistent current 
one finds 
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(22) 



where S s = S s — ha s ir/2. For an ensemble with a large 
dispersion of sizes or kp, only the Cooperon term (s' is 
s time reversed) survives the kp average, and we obtain 



(7(0)) «(/( 2 >(0) 
eh 



,y^ TT^T^ ™fWr)- ( 23 ) 
tttoA ^ |det (M s - I) \ \ (f> J 



By using the H-OdA sum rule and the Gaussian winding 
number distribution, one obtains the average persistent 
current as 



(H4>)) 



V2eh 

ir'l 2 mAu) l 



^/(w c ;n)nsin(47m— ) , (24) 



where the function f(uj c ;n) is defined as 



/(w c ;n) 



dx 5 — exp . 

sinh x \ 2ujf.x 



(25) 



Using the asymptotic expression sinh a; ~ e x /2 valid in 
high-temperature regimes, x = T/T c > 1, this integral 
can be performed analytically. The average persistent 
current for high temperature regime is thus given by 



W)> 



eh 



lAbJr 
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ra=l 



1 H n e 



x sin 47m 



T>T t 
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(26) 



(I (</>)) shows the halving of the flux period with respect 
to the nonaveraged persistent current I{4>), i.e., eq. (14|), 
characteristic of ensemble results. Moreover the ampli- 
tude of (I) for chaotic systems is positive, and scales as 
(I) ~ {kpt) , which is much smaller than that of inte- 
grable systems,E2l (I) ~ kp£. Note that (I(<f>)) decreases 
rapidly with increasing temperature, and the higher har- 
monics die out more quickly. 

In conclusion, we have studied the persistent current in 
quantum chaotic systems by using the Gutzwiller trace 
formula. We derived semiclassical formulas of a typical 
persistent current I typ for single billiards and the aver- 
aged persistent current (I) for an ensemble of billiards 
at finite temperature. Moreover, we have shown that the 
persistent current amplitude scales as P vp ~ kp£ and 
(7) ~ (kp£) . The former result coincides with von Op- 
pen and Riedel's dimensional analysis Hi The different 
kp dependence of persistent current between integrable 
(pyp rv (kpt) ' and (I) ~ kpl) and chaotic systems 
is attributed to the different contributions from noniso- 
lated and from isolated periodic orbits. Therefore, we 
can conclude that the quantum persistent current in bal- 
listic systems becomes a tool to distinguish chaotic and 
regular classical dynamics experimentally. Unfortunately 
no experiment is yet available for chaotic billiards. Ell We 
hope our results will encourage future experiment. 

I would like to thank B. Friedman, Y. Takane, and H. 
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